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APPROXIMATE  ANALYSIS  OF  HEAT  TRANSFER  IN  TRANSPIRED  BOUNDARY 
LAYERS  AT  LIMITING  PRANDTL  NUMBERS 

This  report  presents  a  simple  analytical  procedure  for  approximate 
calculations  of  wall  heat-transfer  rates  in  transpired  laminar 
boundary  layers  at  limiting  Prandtl  numbers.  The  work  represents 
a  further  extension  of  the  investigation  by  Zien  (NOLTR  73-17)  «.n 
heat-transfer  calculations  at  moderate  Prandtl  numbers.  The 
simplicity  and  accuracy  continue  to  be  the  principal  merits  of 
the  method. 
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INTRODUCTION 


A  simple,  approximate  method  for  calculating  the  skin 
friction  on  porous  surfaces  has  recently  been  described  by 
Zien  (Refs.  (1)  and  (2)).  It  is  based  on  an  idea,  due  originally 
to  Volkov  (Ref.  (3)),  for  refining  the  classical  Karman-Pohlhausen 
(K-P)  momentum  integral  technique  in  the  boundary-layer  theory, 
and  is  basically  a  one-parameter  type  of  integral  approach. 

The  chief  merits  of  this  analytical  method,  as  illustrated  in 
References  (1)  and  (2),  lie  in  the  remarkable  combination  of 
simplicity  and  accuracy.  The  relative  insensitivity  of  the  results 
to  the  choice  of  velocity  profiles  is  particularly  noteworthy. 

These  encouraging  findings  are  indicative  of  the  potential  of  the 
method  for  development  into  a  useful  tool  for  practical  calculations 
of  boundary-layer  flows  of  a  more  complex  nature.  Further 
exploitation  and  extension  of  this  new  method  appear  warranted. 

This  paper  presents  a  procedure  for  approximate  calculations 
of  wall  heat-transfer  in  a  transpired  boundary  layer  at  extreme 
Prandtl  numbers.  The  procedure  is  based  on  a  further  modification 
and  extension  of  Volkov's  (Ref.  (3))  original  idea.  Thus,  the 
heat  transfer  on  an  aerodynamic  surface  is  calculated  by  considering 
the  energy  balance  across  the  entire,  attendant  (thermal)  boundary 
layer  at  a  local  flow  station.  A  second  integration  of  the  energy 
equation  in  the  direction  normal  to  the  surface  is  then  performed 
to  provide  an  ordinary  differential  equation  for  the  determination 
of  the  basic  profile-parameter.  To  be  sure,  the  idea  of  using  a 
double-integration  scheme  had  appeared  earlier  in  Whitehead's 
(Ref.  (4))  investigation  of  momentum  boundary  layers.  However,  the 
determination  of  skin  friction  in  Reference  (4)  was  based  on  the 
local  slope  of  the  assumed  velocity  profile,  as  in  the  original 
K-P  procedure. 

For  easy  demonstration  of  the  method,  only  simple,  yet  basic, 
flows  are  considered  in  this  paper.  Thus,  applications  are  made 
to  two-dimensional,  incompressible,  laminar  boundary  layers  of  a 
single-component  fluid.  Results  are  presented  here  only  for 
limiting  cases  of  large  or  small  Prandtl  numbers  where  certain 
simplifying  approximations  can  be  effectively  exploited  to  further 
facilitate  the  exhibition  of  the  central  ideas  of  the  method. 
Detailed  results  for  Pr  =  0(1)  pertaining  to  the  semi-infinite 
flat  plate  with  uniform  surface  mass  flux  are  available  in 
Reference  (5).  The  frictional  heating  is  neglected  in  the  energy 
equation  throughout  the  calculations,  and  this  approximation  is 
generally  valid  for  incompressible  boundary  layers  where  the  Eckert 
number  is  usually  snail.  An  exact  and  comprehensive  analysis  of 
heat  transfer  in  a  class  of  self-similar  transpired  boundary  layers, 
including  the„effects  of  frictional  heating,  can  be  found  in 
Gersten  and  Korner  (Ref.  (6)). 
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The  method  i*  first  applied  to  the  case  of  transpired 
boundary  layers  with  large  Prandtl  numbers.  A  distinguished 
limit  of  large  Prandtl  number  and  small  mass-transfer  rate  is 
readily  apparent  in  the  present  integral  formulation.  Some  examples 
of  this  limiting  flow  are  studied.  The  procedure  is  then  used  to 
study  a  class  cf  low  Prandtl  number  heat-transfer  problems,  and 
the  usefulness  of  the  present  method  in  providing  approximate 
solutions  to  general  transient  heat-conduction  problems  becomes 
eviuent  in  the  process. 

It  is  r.oteci  heic  that  all  the  calculations  presented  in  this 
paper  involve,  at  most,  the  numerical  integration  of  a  single, 
first-order  ordinary  differential  equation.  In  many  cases 
solutions  are  obtained  in  explicit,  analytical  forms.  The  primary 
purpose  of  the  paper  is  not  to  present  solutions  to  any  new  problems. 
Rather,  it  is  to  show  how  simply  some  old  solutions  can  be 
reproduced  with  good  accuracy  by  the  present  method. 

GENERAL  FORMULATION 
Equations  and  3oundary  Conditions 

In  terms  of  the  nondimensional  velocity  (u*,v*)  and  tempera¬ 
ture,  0,  the  basic  differential  equations  and  boundary  conditions 
for  an  incompressible,  constant-property,  laminar  boundary  layer 
over  an  isothermal  porous  surface  are: 
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In  the  usual  practice  of  integral  methods ,  the  asymptotic 
conditions  are  jpecified  at  the  edges  of  the  boundary  layer,  6  and 
6<r»  for  u*  and  0,  respectively.  Some  initial  conditions  on  u*  and 
0  are  generally  required  at  x  ■  0  to  complete  the  formulation. 

These  initial  conditions  are  replaced  by  conditions  on  initial 
boundary- layer  thicknesses,  6(0)  and  6<p(0),  in  the  present  integral 
formulation. 

The  skin-friction  coefficient,  Cf,  and  the  Nusselt  number,  Nx, 
are  related  to  u*  and  0  through  the  following  expressions: 


and 


v  3u* 
u*  ?y“ 


(3a) 

(3b) 


Present  Method  of  Solution 


Since  the  momentum  equation  is  decoupled  from  the  energy 
equation,  the  solution  of  u*  can  proceed  independently  of  the 
temperature  field.  In  order  to  present  a  self-consistent  develop¬ 
ment  of  the  calculation  procedure,  we  shall  make  use  of  the  approxi¬ 
mate  solutions  of  u*  by  a  similar  calculation  scheme  as  reported 
by  References  (1)  and  (2)  for  the  examples  to  be  discussed  below. 

We  now  proceed  directly  to  solve  the  energy  equation.  The 
first  integration  of  (lc)  gives 


k  - 11  -  e!x  /"*d*  ■  t  (If  -  $|  )  (4) 

0  o  w 

where  the  continuity  equation  has  been  used  to  eliminate  v*  in 
favor  of  u*. 

Equation  (4)  is  taken  as  an  expression  for  the  wall  heat  flux 
by  letting  the  integration  cover  the  entire  thermal  boundary  layer. 
Thus  using  (3b)  ,  we  have 

st s  *7?  ■  a;  /T“*e  *y  - c  ,5> 

Equation  (5)  is  a  nondimensional  version  of  the  energy 
balance  across  the  thermal  boundary  layer  at  the  station  x. 


3 


NOLTR  74-121 


'fhcreforo,  the  surface  heat  transfer  is  determined  from  an  integral 
expression  involving  the  assumed  profiles  for  u*  and  0.  Of  course, 
this  representation  would  be  exactly  the  same  as  the  one  using 
the  derivative  of  0  at  the  wall  if  the  profiles  used  for  u*  and  0 
were  the  exact  solutions  of  the  problem  under  consideration. 

The  basic  parameter  of  the  thermal  boundary  layer,  6.r,  is  to 
be  determined  by  the  following  differential  equation  which  results 
from  a  second  integration  of  the  original  energy  equation.  Thus, 
ctfocting  another  integration  of  (4),  we  arrive  at 

4  4  6 

/T  y  f  T  y  j*  T 

dy  f  u*0dy  (c-  ^  /  u*dy)0dy  =  &T  u*0  dy  -  *-  (6) 

In  deriving  (6)  use  has  boon  made  of  the  wall  heat-flux  expression, 
(5)  . 


We  assume  the  following  si.nplc  form  for  6: 

n  i 

0  =  g(nT)  =  £k.n  ,  k.  =  const  (7) 

1  i=0  1  1 

The  choice  of  the  profiles  is  obviously  crude  in  the  sense 
that  they  presume  similarity  for  general  flows  to  be  considered 
below.  However,  they  are  to  be  chosen  such  that  the  essential 
boundary  conditions,  (2b)  and  { 2e ) ,  are  observed. 

Equation  (6)  reduces  to  a  first-order,  nonlinear,  ordinary 
differential  equation  for  6^,  once  a  form  for  0,  such  as  (7), 
is  introduced.  With  6  solved,  the  Stanton  number,  St,  follows 
readily  from  (5)  through  an  algebraic  process. 

APPLICATION  TO  CASES  OF  LARGE  Pr 


General  Remarks 


The  existing  literature  on  boundary-layer  flows  with  large 
Prandtl  numbers  is  voluminous.  Our  purpose  here  is  merely  to 
describe  the  application  of  the  present  simple  method  to  this 
class  of  flows.  Hence,  only  the  references  pertinent  to  the 
present  development  and  results  will  be  mentioned. 

For  large  Prandtl  numbers,  the  simplified  representation  of 
the  velocity  in  the  thermal  boundary  layer  is 


u 


(8) 
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This  approximation  was  introduced  and'  established  by  Fage  and 
Falkner  (Ref.  (7) )  ,  Lighthill  (Ref.  (8)),  Morgan  and  Warner  (Ref. 

(9)),  among  others,  and  is  convenient  for  systematic  developments. 
Subsequent  extensions  and  improvements  of  this  simplifying 
approximation  can  be  found,  for  example,  in  Curie  (Ref.  (10)),  and 
more  recently,  Chao  (Ref.  (ID).  This  simple  expression  will  be 
exploited  in  the  present  formulation.  The  examples  to  be  presented 
here  are  porous-plate  boundary  layers  for  which  the  skin  friction, 
tw,  has  been  calculated  earlier  by  Zien  (Refs.  (1)  and  (2)).  A 
self-consistent  development  of  the  procedure  is  thus  readily 
accessible.  The  results  to  be  presented  and  discussed  in  the 
following  are  all  based  on  the  approximate  solutions  of  tw 
pertaining  to  the  linear  velocity  profile  used  in  Reference  (1) . 
H.vever,  two  temperature  profiles,  g^  and  gj,  are  used  in  the 
calcula'.  ion  to  test  the  sensitivity  of  results  to  the  temperature 
profiles.  We  note  here  that  the  velocity  profile-parameters  take 
the  values  (a,8,y)  -  (1/8,  1/2,  1/6)  in  the  ensuing  calculations 
(see  Ref.  (1)),  and  the  temperature  profile-parameters,  (a,8,v) 
for  g^  and  g^  are  (1/8,  1/2,  1/6)  and  (1/28,  3/10,  1/15),  respectively. 

Porous  Plate:  Similarity  Case 

The  basic  equations  sure  (5)  and  (6)  .  From  (8)  we  have 


u*  -  5  cfRTnT  -  FRTnT 


(9) 


where  R_  ■  u*6_/v.  For  this  particular  flow,  the  result  of  F  as 
reported  by  Zien  (Ref.  (1)),  i.e., 


F  =  Rx"1/2(j  -  BX)’1/2(*  -  D  ,  (10a) 

will  be  used.  Also,  X  is  related  to  X  through 

X  -  X  (£  -  8X)"1/2  .  (10b) 

Substitution  of  (9)  into  (6)  leads  readily  to 


s*r  as;  (BIF>  -  k  *  ®eRt 


(11) 


For  the  case  of  similarity  blowing  (or  suction) ,  eR_  »  const. , 
and  (11)  is  easily  integrated  to  give  the  following  solutions 
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2  ex3 
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lx.-—  1 

I  1  5  Pr 
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<} 


X) 


0,  (12) 


where  X1  =  eR,r. 
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The  heat  transfer  is  obtained  from  (5)  ,  once  R-,  is  found. 
The  result  is 


-  Pr  J  (|)2/3(pp  +  Axe)2/3  (§  -  0X)"1/6  (J  -  X)1/3-  xj  .  (13) 


Equations  (12),  i  1 3 )  and  (10b)  constitute  the  solution  for  N  /R  ly^2 
as  functions  of  X  and  Pr. 

It  is  interesting  to  note  the  behaviof  of  the  present  approxi¬ 
mate  solutions  in  the  limit  of  Pr  ■*  «  and  X  ■  0(1)  for  suction  and 
blowing  separately.  Following  limiting  solutions  are  easily 
deduced  from  (12)  and  (13) : 

11  Nv 

X  <  0  :  Xx  -  -  i  ^  and  — ^  -  |  X  | Pr  ;  (14) 

e  Rx 


X  >  0 


X'3/2 


-  OX) 


1/2' 


nl/2 


i- 


and 


— U  -  -  pr  X  (I  B  -  1)  .  (15) 

■ST 

For  the  special  case  of  X  -  0,  the  asymptotic  limits  for 

Nx/Rx*^2  as  Pr  -*■  •  are  0.367  Pr*^3  and  0.338  Pr*^3  for  0  »  g^  and 

g.,  respectively.  These  approximate  limits  are  to  be  compared  with 
’  1/3 

the  exact  linut  of  0.339  Pr  (see,  for  example,  Ref.  (12)). 


As  Pr  ■*  «,  the  heat-transfer  rate  predicted  here  approaches 
the  correct  asymptote  (see  Gersten  and  Korner,  Ref.  (6))  independent 
of  profiles  for  the  case  of  suction.  For  blowing,  the  present 
solutions  show  a  strong  profile-dependent  behavior.  Note  that  the 
quantity  in  the  parenthesis  of  (15)  is  negative  for  both  profiles 
considered  here.  The  implications  of  the  above  observations  are 
thus  clear.  As  Pr  increases,  the  present  solutions  are  uniformly 
good  for  suction,  but  predict  too  capid  a  decay  in  the  heat-transfer 
rate  for  blowing,  especially  when  X  is  large. 
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Some  typical  results  are  illustrated  in  Figures  1-4  and 
compared  with  the  available  exact  solutions  due  to  Stewart  and 
Frober  (Ref.  (13) ),  and  Thompson  (Ref.  (14))  whenever  appropriate. 

Figure  1  shows  the  results  for  Pr  ■  1  for  which  on  exact 
integral,  9  ■  1  -  u*,  holds.  Exact  solutions  for  heat  transfer  are 
thus  available  from  the  exact  skin-friction  results  of  Emmons  and 
Leigh  (Ref.  (15)),  or  directly  from  Reference  (13),  for  comparison. 
It  is  remarkable  that  the  present  procedure,  suitable  for  large 
Prandtl  numbers,  yields  very  accurate  results  of  heat  transfer 
for  Pr  -  1.  Of  course,  the  finding  that  the  validity  of  a 
calculation  procedure  designed  for  large  Prandtl  numbers  actually 
extends  into  regions  of  moderate  Prandtl  numbers  is  not  totally 
unexpected,  (see  Ref.  (8)).  We  also  note  the  insensitivity  of 
the  present  results  to  the  profiles  used  in  the  calculations. 

Figures  2-4  show  the  heat-transfer  results  are  functions  of 
the  Prandtl  number  at  three  values  of  X:  and  -  5, 

/I  /2  10 

for  which  exact  solutions  are  available  for  comparison.  These 
results  clearly  confirm  our  aforementioned  predictions  based  on 
the  ~ehavior  of  the  present  solutions  in  the  limit  of  Pr  -►  <■». 
However,  it  should  be  remarked  here  that  aa  Pr  ♦  ■,  the  results 
pertaining  to  blowing,  though  quantitatively  inaccurate,  are 
still  qualitatively  correct.  It  is  possible  that  the  thermal 
boundary  layer  is  blown  off  the  surface  while  the  momentum  boundary 
layer  still  remains  attached. 

The  same  approach  and  profiles  are  then  employed  in  the 
original  K-P  procedure  for  this  problem,  using  tw  obtained  from 
the  original  K-P  method  with  a  linear  velocity  profile,  as 
reported  in  Reference  (1).  The  results  are  presented  here  only  for 

A  »  — —  in  Figure  5.  A  comparison  between  Figure  5  and  Figure  3 

/? 

clearly  illustrates  the  superiority  of  the  refined  K-P  procedure. 
The  absence  of  a  sharp  decrease  of  the  heat-transfer  rate  in  the 
interval,  6  <  Pr  <  100,  is  typical  of  the  failure  of  the  original 
K-P  procedure.  Nevertheless,  it  should  be  mentioned  that  for 
the  case  of  suction,  the  original  K-P  procedure  with  the  same 
profiles  is  also  capable  of  predicting  the  correct  asymptotic 
behavio”  for  Pr  ♦  »,  only  the  sensitivity  of  results  to  profiles 
is(  more  pronounced. 

A  Distinguished  Limit 

The  singular  perturbation  nature  of  the  boundary-layer  energy 
equation  in  the  limit  of  Pr  +  »  is  well-known.  (See,  for  example, 
Narasimha  and  Vasantha,  Reference  (16).  It  is  also  apparent  in  (lc) 
in  which  the  coefficient  of  the  highest  derivative  approaches  zero 
as  Pr  +  ».  The  asymptotic  structure  of  the  thermal  boundary  layer 
with  surface  blowing  in  the  limit  of  Pr  *  has  been  discussed 
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by  Kassoy  (Ref.  (17)),  among  others.  The  interesting  phenomenon 
of  thermal  layer  blowoff  is  qualitatively  analyzed  for  the  case 
of  a  porous  plate  with  similarity  blowing  in  Reference  (17). 

Physically  speaking,  the  thermal  boundary  layer  has  a 
vanishing  thickness  relative  to  the  momentum  boundary  layer  as 
Pr  -*  ®,  because  the  viscous  diffusion  (~v)  dominates  the  thermal 
diffusion  ('<) .  Therefore,  the  heat  transfer  takes  place  in  a 
very  thin  region  close  to  the  wall  where  the  convective  velocity 
components  (u,v)  are  also  small.  In  the  case  of  transpired 
boundary  layers,  an  additional  mechanism  of  heat  transfer  (~e) 
appears.  The  precise  orders  of  magr itude  of  these  physical 
quantities  as  Pr  ■*  ®  are  related  uniquely  only  in  a  distinguished 
limit  in  which  convection,  conduction  and  mass  transfer  are  of 
equal  importance.  It  is  convenient  and  instructive  to  deduce 
such  a  limit  from  the  equations  derived  in  the  process  of  the 
present  development.  We  remark  here  that,  as  Pr  +  «°,  the  present 
method. seems  particularly  suited  for  this  case  of  vanishingly 
small  A  in  view  of  the  results  discussed  in  the  section  of  Porous 
Plate:  Similarity  Case. 

Let  us  consider  (11) .  The  three  terms  from  left  to  right 
represent,  respectively,  the  effects  of  convection,  conduction 
and  mass  transfer.  Therefore,  we  require  that,  as  Pr  -*•  °°, 


i_ 

Pr 


eR„ 


(16) 


-1/2 

where  we  have  used  the  relation  F  *  0(Rx  ). 


It  follows  immediately  that 


Pr  •*  ®:  A 


0(Pr_2/3) 


'  X1 


0  (Pr-1)  and  6^/6 


0(Pr"1/3) . (17) 


Thus,  the  distinguished  limit  corresponds  to  vanishingly 
small  values  of  A  and  6_/5.  We  are  now  led  to  introduce  the 
following  variables  forthe  study  of  this  limiting  flow: 


A  «  Pr2/3  A  =  0(1) 

(18a) 

Ax-  Pr  Ax  =  0(1) 

(18b) 

In  this  limit,  F  may  be  expanded  as 

F  =  Fg  +  A  F^  +  . . . 

(19) 
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where  Fg  corresponds  to  the  skin-friction  coefficient  in  the 
absence  of  mass  transfer.  Only  the  leading  term,  Fg,  will  be 
needed  in  the  ensuing  calculations.  The  result  to  be  obtained 
thus  corresponds  to  an  approximation  of  the  leading  asymptotic 
term  in  this  distinguished  limit.  A  systematic  asymptotic  analysis 
of  the  general  flow  in  this  distinguished  limit  appears  to  be  an 
interesting  problem  of  fundamental  importance  but  beyond  the  scope 
of  the  present  paper. 

For  a  flat  plate,  Reference  (1)  provides 


and,  as  mentioned  earlier,  (a,y)  ■  (1/8,  1/6)  is  used  in  the 
calculations. 

Porous  Plate:  Similarity  Case 

Introducing  A  and  Aj,  and  applying  the  limit  process,  (17),  to 
(12)  and  (13),  we  arrive  at  the  following  results: 

A3  -  (2)3/2(c)1/2A3(1  +  AjBJ/lay)  (21) 

and 

-  Pr1/3  [y  (a)“1/6  (5)  "2/3  <y)  1/3  (2)“1/2  (1  +  A10)2/3  -  Aj  .(22) 

Results  of  Nx/(Rx1//2Pr1//3)  as  a  function  of  A  for  the  two 

temperature  profiles  are  plotted  in  Figure  6.  It  can  be  shown 
that  both  curves  approach  the  exact  limit  of  asymptotic  suction. 

The  dependence  of  the  results  to  profiles  is  still  weak  in  this 
limiting  flow,  and  the  thermal  boundary- layer  blowoff  (nx  *  0) 
is  predicted  at  a  finite  value  of  A.  An  exact  treatment  of  this 
limiting  case  was  recently  reported  by  Gersten  (Ref.  (18))  for  a 
class  of  Falkner'-Skan  flows.  The  same  distinguished  limit  was 
derived  in  a  different  manner  in  Reference  (18) .  The  exact  solution 
pertaining  to  the  flat-plate  configuration  is  included  in  Figure  6 
for  comparison  with  the  present  approximate  results.  The  agreement, 
in  general,  is  clearly  quite  satisfactory,  except  near  thermal 
boundary-layer  blowoff  where  the  exact  solution  exhibits  an 
exponential  decay  in  heat  transfer. 

Porous  Plate:  Nonsimilar  Case 


Equation  (11)  is  still  the  basic  differential  equation  for 
this  nonsimilar  flow,  as  the  velocity,  u,  is  approximated  by  (8). 
Now,  applying  the  limit  process,  (17)  and  (18),  to  (11)  and  noting 
(20) ,  we  have  the  following  differential  equation  for  the  limiting 
flow: 
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and 


dAj 

a r 


(8a) 


1/2 


.  _1£L 


(1  + 


SA1) 


A1 

<r> 


h2 

2<7T> 


A 


i8J*l 

(ya) 


1/6 

T7I 


(23a) 


(23b) 


Substituting  (9)  into 
obtain 


(5)  and  introducing  A  and  A^,  we 


(24) 

Results  of  N  / (R^^^2Pr^^  3)  as  a  function  of  A  are  plotted 

in  Figure  7  for  the  two  temperature  profiles  used  in  the 
calculation.  The  behavior  of  heat  transfer  is  cuite  analogous 
to  that  for  the  similarity  case,  Figure  6.  We  also  note  that 
the  two  curves  will  converge  in  the  region  of  large  negative  A 
and  approach  the  exact  asymptotic  suction  limit, 

N  /(R  1/2pr1/3)  -  | A | .  No  other  limiting  solutions  to  this  non¬ 
similar  ilow  seem  to  exist  in  the  literature,  to  the  best  of  the 
author's  knowledge. 

APPLICATION  TO  CASES  OF  SMALL  Pr 


172 


Pr 


1/J 


i 

a 


V  “ 


-  !*)a 


General  Remarks 


In  this  section,  the  method  is  applied  to  the  calculation 
of  boundary- layer  heat  transfer  in  the  limit  of  Pr  ♦  0.  The 
mathematical  limit  here  is  R  ■*  ®,  Pr  •+  0  and  R  Pr  -*■  »,  so  that 
the  concept  of  a  thin  thermal  layer  near  the  surface  is  still 
valid.  The  chief  simplification  is  that  the  velocity  profile,  u, 
inside  the  thermal  boundary  layer  can  be  approximated  by  the 
external  potential  stream  ue (x) ,  (Ref.  (12)).  Therefore,  we  need 
only  to  assume  the  temperature  profile.  The  resulting  solution 
corresponds  to  the  limiting  case  of  Pr  -  0,  i.e.,  an  inviscid 
but  heat-conducting  fluid.  Three  representative  examples  will  be 
presented  in  the  following.  These  examples  will  explicitly  reveal 
the  relevance  of  the  differential  equations  to  a  class  of  heat- 
conduction  problems.  Therefore,  the  success  of  the  present  method 
here  would  imply  its  potential  usefulness  in  studying  general 
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problems  in  transient  heat  conduction  with  complex  boundary 
conditions,  including  nonlinear  cases.  Some  initial  attempts 
have  already  been  made  by  Volkov  and  Li-Orlov  (Ref.  (19)).  A 
more  thorough  investigation  of  this  application  appears  desirable, 
and  should  lead  to  a  modified  version  of  the  work  of  Goodman 
(Ref.  (20))  who  had  applied  the  original  K-P  method  to  a  variety 
of  heat-conduction  problems. 

Porous  Plate:  e  ■  const.  -  0 (l//RxPr) 

Here  the  velocity  (u,v)  -  (u«,vw)  and  the  energy  equation 
takes  the  form 

55  *  E  3y  ’  17  •  ,25) 

The  first  and  second  integrations  of  (25)  lead  to,  respectively, 

Nx  -  dl>T 

V*  ■  “i  as;  -  c  1261 


dR 

li  as1- 


PrRT 


+  eB 


where  (d.,B)  aie  profile-parameters  defined  in  the  Nomenclature. 

In  deriving  (26)  and  (27)  ,  a  temperature  profile  of  the  form 
0  *  g (Ht^  is  assumed.  Equation  (27)  is  readily  integrated  to  give 


a.  a. 

c2RxPr  -  ^  ER^r  -  In  (1  *  BeRTPr) 
and  the  heat-transfer  coefficient  follows  from  (26)  as 


/TTr  “1 


(1  -  ^-)e/PrR" 


Equations  (28)  and  (29)  form  a  parametric  representation  for  the 
heat-transfer  coefficient,  i.e.. 


f  (e/FPPr ) 
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Results  corresponding  to  the  two  temperature  profiles, 
and  g a,  are  illustrated  in  Figure  8.  The  profile-parameters 
associated  with  gj  and  g4  and  (a,,f5)  -  (1/3,  1/2)  and  (2/15,  3/10), 
respectively. 

It  is  interesting  to  note  that  (25)  is  also  relevant  to  the 
problem  of  the  transient  heat  conduction  in  a  semi-infinite  solid 
moving  with  a  constant  velocity  c  and  maintained  at  a  constant 
temperature  at  y  -  0.  The  exact  solution  is  available  in  Carslaw 
and  Jaeger  (Ref.  (21)).  The  heat-transfer  rate  is  readily  obtained 
from  the  solution  for  0,  i.e.. 


Nv  ,  c2R  Pr  c. /R  Pr  c /R  "Pr 

- X—  *  —  exp ( - T^—) - y—  erfc  (  ■  *X~  )  .  (31) 

✓  RxPr  /n  * 

The  exact  solution  is  also  shown  in  Figure  8  for  comparison 
with  the  approximate  solutions,  (30).  The  present  results  are 
clearly  very  satisfactory  in  accuracy;  the  close  agreement  between 
the  results  corresponding  to  the  two  temperature  Drofiles  is  also 
evident.  However,  it  should  be  pointed  out  that  as  the  blowing 
intensity  increases,  the  present  approximate  solutions  predict 
zero  heat  transfer  (thermal  boundary-layer  blowoff)  at  finite 
values  of  c/R  Pr ,  while  the  exact  solution  shows  exponential  decay 
in  heat  transfer,  lherefore,  as  is  usually  the  case,  the  accuracy 
of  the  present  solution  begins  to  deteriorate  as  the  blow-off 
point  is  approached. 

Circular  Cylinder 

We  now  consider  the  case  of  a  circular  cylinder  in  crossflow. 
The  orthogonal  curvilinear  coordinate  system,  (x,y),  has  its 
origin  at  the  forward  stagnation  point  with  x  measuring  the 
distance  along  the  surface  (x  =  a$).  Then  u*  =  2sin$,  v*  *  -2y* 
cos4».  The  energy  equation  and  the  boundary  conditions  are 

,99  *  4  38  i  a2e  . 

sin*  7*  -  y*  cos»  37*  -  i^pF  —  (32a) 

0(*,  0)  =  1  0  ($,  5*)  =  0  ,  (32b) 

* 

where  u*  »  u/u,,,,  v*  =  v/um,  y*  =  y/a  and  6T  =  6T/a. 

Applying  the  present  method  to  (32a)  results  in  the  following 
two  equations,  based  on  the  assumed  profile,  0  =  g(nT): 
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1  ND  '  d6T 

2  R^Pr  "  2(sln*  W  +  6T  cos*) 

and 

*  d6T  *2  i  i 

5T  sin*  3T  +  4T  cos*  *  57  R^Pr 

with 

«J  (0,o)  *  (OjR^r)-1/2  =  (6*> 


Equation  (34a)  can  be  easily  integrated  to  give 


(33) 


(34a) 


(34b) 


Equation  (33)  then  gives 


16*)  sec  l  .  (35) 

1  0  1 

2  -L—  cos  4  (36) 

✓di  2 

1.73  cos  |  for  6  =  gx  (37a) 

1.64  cos  ^  for  0  =  94  •  (37b) 


The  exact  solution  to  this  problem  can  be  found  in  Grosh 
and  Cesi  (Ref.  (22))  who  transformed  the  equation  to  a  standard 
form  of  one-dimensional  heat  equation.  The  heat-transfer  result  is: 


•V 


=  1.60  cos  j 


(38) 


The  present  method  is  again  shown  to  yield  satisfactory  results. 


Two-Dimensional  Stagnation  Point,  Uniform  Vw 


As  a  last  example  of  this  limiting  flow,  we  consider  the 
heat  transfer  near  a  plane  stagnation  point  in  the  presence  of 
uniform  surface  mass  flux.  Of  course,  the  treatment  is  still 
within  the  framework  of  thermal  boundary- layer  theory. 
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Here  we  have  u 


equation  and  the  boundary  conditions  are 


and  v  ■  Vw  -  yuj.  The  energy 


*  ao  A  ~  ao  aJe 

x*  57*+  <‘l  —*> 


0  (x*,0)  =  1,  0  (x*  ,  ®)  =*  0, 


(39a) 

(39b) 


where  the  coordinates  and  velocity  are  normalized  by  the 
characteristic  length  /k/u,  and  velocity  /tcu,  respectively. 
In  particular,  we  have 


w 


'KUl 


V  /  u.x 

^V—  - c  o(i) 


(40) 


Using  0  =  g  (n<p)  and  applying  the  present  method  to  equation 
(39a)  result  in 


d  6  w 

51**S?  *  ai6T!  ■  1  +  V  {T 


(41) 


and 


Nx 


% 


d6* 

—  -  e  («$  +  x*  xjr- )  -  h  (42) 


/RxPr  K  ‘  W1*  T  - 

The  solution  to  equation  (41)  is  easily  obtained  as 


6* 

T 


±  ['  :i  +  ] 


*  const. 


(43) 


Equation  (42)  then  gives 


— -  =  0  6*  -  A.  (44) 

/RxPr 

The  exact  solution  to  equation  (39)  is  also  easily  obtained  by 
noting  that  the  boundary  conditions  and  the  original  partial 
differential  equation,  (39a)  ,  admit  a  simple  solution  of  the  form 
6=6  (y*) .  Hence,  (39a)  is  essentially  an  ordinary  differential 
equation  whose  solution  has  a  simple  closed-form.  The  exact  heat 
transfer  follows  from  the  Solution  for  0,  and  is  given  by 
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r 


(45) 


The  approximate  solutions  have  been  obtained  for  two  temperature 
profiles,  g,  and  g,,  and  compared  with  the  exact  solution, 
equation  (45),  in  Figure  9.  Again,  the  approximate  solutions  are 
shown  to  be  very  accurate.  It  can  further  be  shown  analytically 
that  the  approximate  solutions  reduce  to  the  exact  asymptotic 

limit  as  •*  -«,  independent  of  profiles.  Near  thermal  boundary- 

layer  blowoff,  the  comparison  resembles  closely  that  of  the  case 
of  a  porous  plate  studied  in  the  Porous  Plate  section  of  this 
report. 


CONCLUDING  REMARKS 

The  primary  purpose  of  this  paper  is  to  present  a  simple 
procedure  for  practical  calculations  of  heat  transfer  in  boundary 
layers.  Oversimplified  and,  in  fact,  improper  profiles  are 
deliberately  used  in  the  calculations,  and  the  principal  merit  of 
the  method,  namely,  the  remarkable  combination  of  simplicity  and 
accuracy,  is  amply  demonstrated  and  fully  exploited.  No  effort  is 
expended  here  in  providing  the  details  of  the  flow  field. 

Results  of  these  calculations  indicate  that  the  rudimentary 
application  of  the  method,  as  described  in  this  paper,  is  adequate 
for  most  engineering  purposes.  The  weak  dependence  of  the  results 
to  the  profiles  continues  to  prevail.  This  is  believed  to  be  a 
consequence  of  using  an  integral  expression  for  the  surface  heat 
flux;  the  effect  of  improper  profiles  is  mostly  reflected  in  the 
somewhat  spurious  predictions  of  transverse  scales  of  the  boundary 
layers,  6  and  6™,  so  that  the  predictions  of  the  wall  heat-flux 
remain  reasonably  accurate.  However,  near  the  thermal  boundary- 
layer  blowoff  where  the  heat-transfer  rates  are  small,  the 
results  are  useful  only  for  qualitative  purposes.  This  finding 
is  reminiscent  of  the  difficulty  encountered  in  the  previous  skin- 
friction  calculations  (Refs.  (1)  and  (2)). 

The  possibility  of  using  the  method  to  study  general  problems 
of  transient  heat  conduction  appears  promising  and  warrants 
further  investigation.  In  particular,  the  present  mi  thod  seems 
capable  of  providing  simple  approximate  solutions  to  heat  conduction 
problems  with  phase  changes  and  moving  boundaries,  such  as  melting, 
ablations,  etc.  These  problems  have  the  characteristic  feature  of 
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nonlinear  boundary  conditions  imposed  on  a  (unknown)  moving 
boundary.  The  technological  importance  of  such  problems  is  self- 
evident.  The  present  approximate  method  should  serve  *s  a  useful 
guide  to  the  development  of  extensive  numerical  efforts  in  this 
area.  Results  of  this  study  will  appear  in  a  separate  report. 
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FIG.  1  POROUS  PLATE,  SIMILARITY  CASE,  Pr  *  1  (PRESENT  METHOD,  u  *yr J\i) 
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FIG.  4  POROUS  PLATE,  SIMILARITY  SUCTION  (PRESENT  METHOD,  u  •  yr  /p> 
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FIG.  8  POROUS  PLATE,  UNIFORM €,  Pr-H) (PRESENT  METHOD) 
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FIG.  9  PLANE  STAGNATION  POINT.  UNIFORM  Vw.  Pr— 0  (PRESENT  METHOD) 


